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Abstract 

The ground-state energy, the effective mass and the number of virtual phonons of 
the optical large polaron confined strictly in one dimension have been estimated by 
using the generalized Gaussian approximation. The leading-order terms take care of 
all Gaussian fluctuations in the system and improve the conventional variational esti- 
mates at finite coupling. Particularly, the lowest upper bound to the polaron ground- 
state energy has been obtained. The non-Gaussian contributions systematically correct 
the leading-order approximations. We have obtained exact analytical solutions in the 
weak- and strong-coupling limit and reasonable numerical data for intermediate cou- 
pling. Our result for the number of excited phonons limits the validity region of the 
few-phonon approximation methods. 



PACS number(s): 31.15.Kb; 63.20.Kr; 71.38.+i; 73.20.Dx 



1 Introduction 

A conducting electron in a polar crystal polarizes the surrounding lattice and, by interacting 
with the quanta of lattice vibration, phonons, forms a quasi-particle, the polaron. 

The polaron problem traditionally has been investigated in three dimensions jl], @j, al- 
though the electron motions may be localized in reduced dimensionality || . Many theoretical 
investigations have been made to study the polaronic properties in lower dimensions (see, 
e.g. 0, U). A considerable research attention has been devoted to the study of the polaron 
confined in one dimension @, 0, 0- In particular, recent studies on the electron-phonon 
interaction in quasi one- dimensional systems have explored a polaronic effect pronounced 



stronger than in two-dimensional structures [10]. There are two distinct models employed 



for the one-dimensional polaron: the small and the large polaron. Small polarons are more 



strictly bounded and observed in linear conjugative organic polymer conductors |TT[] while 
the latter model describes electrons confined in one-direction in polar semiconductors JT2 . 
The phonons have been considered in the optical and acoustic modes, depending on the 
specific dispersion rule 0. Regardless of the extremely confined structure of the system, 
polar optical phonons are treated in the spirit of bulk theory. 

The resulting interaction of an electron confined in one dimension with the phonons of 
longitudinal optical lattice vibrations can be modelled by the following Hamiltonian || [7], || 

H=^- + huj Y,aW + gY,( a k eikr + a U~ ikr ), 9 = huj (ha/muJoL) 1/2 , (1) 
2m V V v ' 
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where m, p and r denote the electron bare mass, the momentum and position of the electron, 
L is the length of lattice crystal, whereas k, and a\ are wave vector, annihilation and 
creation operators of a phonon. The optical polaron model imposes the phonon frequency 
ojq independent on k. Further we will use appropriate units, such that, m = u = h = 1. 
The electron-phonon coupling factor g behaves g oc ^Ja/Q/\k\ ( - d ~ 1 ^ 2 in d dimensions |T3"|, |T4]|, 
where Q is the volume of the system. Therefore, g does not depend on k in one dimension 
[|. The Frohlich dimensionless coupling constant a takes values ranging from 0.02 to 4 
for most actual crystals. In this paper the coupling constant a is identical to 27ra op in ||, 
a' in J7| and a in g. 

Many approaches have been developed to investigate the quasi-particle properties of the 
system ([!]) in the ground state || [F|. Strict results for these quantities are available only in 
the weak- (a — > 0) and strong-coupling (a — > oo) limit (see, e.g. @, ||). 

In the present paper we have made a systematic investigations of the large optical polaron 
confined strictly in one dimension within a path-integral (PI) method. We have obtained 
exact analytical solutions of the self-energy, the effective mass and the number of virtual 
phonons in the weak- and strong-coupling limit as well as reasonable numerical data im- 
proving known results || |7|, § at finite a. Particularly, our result for the number of excited 
phonons limits the validity region of the few-phonon approximations || and shows that as 
a grows, the efficiency of the general Gaussian approximation drops faster than in higher 
dimensions. 



2 Polaron Ground-state Properties 

1. The phonon variables in (HD may be exactly eliminated by using the Feynman PI technique 
0. The partition function for a resting polaron reads 

r(/3)=0 

Z,(a)=e-^ = (e^) Q , <(...)> = JSre-^%..), (2) 

r(0)=0 

1 ? d 2 
S [r] = -J dtds r(t) Do\t, s) r(s) , D \t, s) = -5{t - s) — , 

o 

where the inverse temperature f3 = 1/kgT is infinite in the ground state: f3 — > oo. Therefore, 
E{a) = F(a, oo) represents the polaron ground-state energy (GSE). The polaron retarded 
self-interaction is given by 

P oo 

W[r] = -j= jjdtds e-l'-*! J g e* kR ^ = Jd{l tsk e lkR ^ , R(t, s) = r(t) - r(s) . (3) 

V -oo 

2. For a slow polaron the partition function projected at small fixed momentum p is 

p(f3)=x 

Zb{oc,p) = e ~ ma ' p) = [dxe~ ipx hpe' So[p]+aW[p] . (4) 

p(0)=0 

Since the polaron action is translationally invariant, the energy E(a,p) is a continuous 
function of p. By expanding it around small momentum we define the effective mass (EM) 
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of the polaron m*{a) as follows: 

2 

E ^P) = E (®) + ^a)+0(P 4 )- (5) 

With the new integration variable r(t) = p(t) — q(t), where q(t) = xt//3 = yt, one goes to 
conventional closed-end paths r(t) and rewrites (§]) as follows: 

pE(a,p) _ J J±_ lj„ l(} -iPpy-(Jy 2 /2 -(3T{a,y 2 ) ~^(a,y 2 ) _ / „aW[r+q] 




For (3 — > oo the integral over y in @ may be taken by the saddle point method. The saddle 
point y = satisfies the equation 

ip + t 1 + 2 0] = , ^(«, = ^(a, . (7) 



Then, 

E{a, p) = ip 41 + \ + ^(a, = ^(«) + 2(1+ f^ (a)) + 0(p 4 ) , (8) 
where J-(a) = J-(a, 0) and ^(a) = JF^(«,0). Comparing fl5|) with (J|) we find 

£(a) = .F(a), m*(a) = 1 + 2^(a) . (9) 
To evaluate ^(a) we use the expansion 

W[r + g] = W[r] + VL[r] - y K 2 [r] + 0(y 3 ) , 

I4[r] = jdn tsk e lkR ^k(t-s), V 2 [r] = Jdn tsk e lkR ^k 2 {t-sf . 
Neglecting terms ~ 0(y 3 ) we obtain 

e H«W) = Z/j(a) exp(^|g} , Y,(a) = \{e^ (V?[r\ + V 2 [r])) Q . (10) 



Hence, 



3. The average number of virtual phonons (ANVP) excited by the electron-lattice inter- 
action is given by 



#(a) = <0| 2>to*|0) 

k 



dH 



\ 



. (5 du { 

U) =l ' 



1 °' lnTre-^ 



,9^0/ 

The differentiation in (|I^) should not affect the interaction term in H. We obtain 



(12) 

0)0 = 1 



N(a) = drre-T dk ^ff^ , X,(a,k,r) = ( e ^(M + ^M\ , (13) 



where r = t — s. Note, the following relation between the ANVP and GSE takes place [|T^] 

3 d_ 

2 a d^ 



N{a)=(l-la4-)E(a). (14) 
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3 Generalized Gaussian Approximation 

For finite a, the nonlocality and (^-singularity arising in (^) prevent any explicit evalua- 
tion of Z ^(a), Yp(a) and Xp(a,k,r). Below we represent a nonvariational approximation 
method, which isolates completely the Gaussian contributions from these Pis. The remain- 
ing non-Gaussian corrections turn to be relatively small (at least for the energy) and may 
be systematically estimated to improve the leading-order approximations. 

i. First, we demonstrate the basic idea of our method by evaluating Z 13(a). Note, the 
initial representation @ is optimal only for a — > 0. To describe the system more efficiently 
at finite a we introduce a functional integration weighted with the most general Gaussian 
measure da[r] as follows 

r(/3)=0 

((...))= J da[r] (...), da[r} = 5re-^ r ' D ~ lr \ (1) = 1 . (15) 

r(0)=0 

We assume that all Gaussian configurations of the system is totally concentrated in Ax[r]. 
This imposes specific requirements on D~ l (t, s) (and on its Green function D(t, s)) that will 
be discussed later. Obviously, the following relations take place 

(r(t) r(s)) = D(t, s) , (e ikR ^) = e~* F <^'V , F(t, s) = + D ^^ S ^ _ D ( t) s ) . 



Then, we rewrite @) in the new representation as follows 



( 1 i) 1 / 11 

Z p (a) 



exp < -Tr In — > ( exp <^ - // dtds r(t) D~ l (t, s) r(s) - S [r] + a W[r] 



(16) 



Let us introduce the conception of the normal- ordered form of functionals with respect 
to dcr[r]. In particular, we use the following normal forms: 

:r (t) r(s):=r(t) r(s) -D(t,s), : e ikR{t ' s) - = e ikR{t > s) e k2p{t ' s) (17) 

so that 

(:r(f)r(a):> = 0, ( :e ^ (M) :) = 0, e\ = e* - 1 - z - y . (18) 
Functional W[r] may be decomposed as follows: 

W[r] = JdA tsk :e tkR ^ : = W[0] + iW l [r] - ^W 2 [r] + W int [r] , 

W 1 [r} = JdA tsk k:R(t,s):, W 2 [r] = JdA tsk k 2 :R 2 (t, s) : , (19) 

W int [r] = JdA tsk : ef R{t ' s) : , dA tsk = dQ tsk e ~ k2p ^ . 

Since all quadratic configurations of the polaron action in the new representation is totally 
included in da[r], any extra quadratic parts should be eliminated as follows 

-JJdtds:r(t)[D~ 1 (t,s)-D 1 (t,s)]r(s):+aW 2 [r} = 0, Vr . (20) 
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This requirement leads to the following constraint equations for function F(t): 



F(t) 
D(k) 



7T 



- dk [1 - cos(fct)] D(k) , 





k 2 



a 



/27T 



dt e 



1 — cos(fct) 
F 3 / 2 (t) 



(21) 



where D(k) is the Fourier transform of D(t). Note, Do(k) = 1/k 2 . Taking into account (|20|) 
we rewrite (|I6|) 



Z /9 (ot) 
E (a) 



1 

~2^ 

o 



In (A; 2 D(A;) 



^(ot) = (e 
fc 2 D(A;) + 1 



aW int [r] 



a 



/27T 



exp(— £) 
FV2(t) 



(22) 



where i?o(ai) is the Gaussian leading-order contribution to the GSE of the polaron confined 
in one dimension. The non-Gaussian corrections associated with Jp(ot) can be estimated 
systematically. Equations (|T6"1) -(E22"|) serve as the basis of the Generalized Gaussian Approxi- 
mation (GGA) method applied to Zp(a). 

Note, E (a) represents a upper bound to the true GSE of the one-dimensional polaron. 
Indeed, (Wj n i[r]) = 0- Then, by using the Jensen-Peierls inequality one obtains 



M<*) > e a{w '" t[r]) 



so, 



E(a) < EM . 



(23) 



ii. Now we evaluate Yp(a) within the GGA method. Remember, transformations (K)- 
([22] ) affect only the exponent in (pT0|) , but not any prefactor of the exponential. Therefore, 
despite the factor (Vj 2 [r] + V 2 [r]) /2, all steps of the GGA method may be directly applied 
to Yp(a). Constraint equations ( pl|) remain without any change. We obtain 



where 



-f3Eo(a) _ ^ / aW int [r] 



(V?[r] + V 2 [r])) = e-^W faa) + \m<*)} 



IfdAtsk k 2 (t-s) 
dA tsk k 2 (t-s) 2 ( <> ik ; 



J /3 (a) = (e aW/ -M-Q[r]) , 



(24) 



(25) 



1 + 



rfA tsfc A;(t- S )(e ifc ^ 



Q[r] 



Substituting (p4|) into (JXTJ) and @ we obtain the leading-order (Gaussian) approximation 
to the polaron EM as follows 



mS(a) = l+2^(a) = l + -^=yi 

o 



F 3 / 2 (t) ' 



(26) 



Higher-order (non- Gaussian) corrections to the EM are given by 

Ma) 



Am*(a) 



lim 



2^ ) (a)-[J /3 (a)-l] . 



(27) 
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iii. Note, the additional linear term ikR(t, s) standing in the exponent in (|T3|) affects 
neither the normal ordering, nor the elimination of the quadratic path configurations. Hence, 
the constraint equations ( f2~T|) do not change for the ratio Xp(a, k,t)/Zp(a). So, we write 



X/3(a ' fc,r) = e~ k2p ^ + non-Gaussian part . (28) 

Substituting (p8|) into flllf ) and neglecting the non-Gaussian part we obtain the leading- 
order approximation to the ANVP as follows 

oo 



No{a) = v^J dt ^m- (29) 







4 Non-Gaussian Corrections 

To improve the obtained Gaussian approximations we use the following expansions 

oo n oo n 

M*) = E ~J W5M , M«) = E ZJ (Wnt[r] QM) ■ (30) 



n=2 711 "0 U - 



These are not plain perturbation expansions over a because each of terms in angle brackets 
contains a in a complicated way so that these series converge fast even for large a. Taking 
into account higher-order non-Gaussian contributions, we obtain corrected values 

E n {a) = E {a) + AE 2 {a) + . . . + AE n {a) , 

m* n {a) = m* (a) + Am* 2 (a) + . . . + Am* n (a) , (31) 
N n {a) = N {a) + AN 2 {a) + . . . + AN n {a) . 



In particular, 



oP'ti ffff e I* s ^ \ x ^ 

AE2ia) = -WJJf tdSdXdy p^s,x,y) ' 

oP"j[ rrrr g — I* - s \~ \ x ~ y\ r 

Am* 2 (a) = — JJJJdtdsdxdy p3/2{t s ^ ^ {{t - s) F(x - y) 

+ {x-y) 2 F{t-s) + \t-s\\x-y\E{t,s,x,y)} , (32) 

p(t,s,x,y) = 4 F{t - s)F(x - y) - E 2 (t, s,x,y) , 
E(t,s,x,y) = F(t-x) + F{s-y)-F(s-x)-F(t-y). 



Note, the integrations over s and y in (|32D may be taken explicitly. The remaining double 
integrals can be derived analytically in the weak- and strong-coupling limit, and numerically 
for finite a. 

5 Exact and Numerical Results 

We admit that the accuracy reached with the second- and third-order non-Gaussian correc- 



tions in (31) is sufficient for our consideration. Below we have represented exact analytical 



solutions in the weak- and strong-coupling limit and reasonable numerical data for finite a. 
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Weak coupling 

For a — > the solution of the constraint equations fl2~i|) reads 



D(k) = \k 2 + 2a(V2\jl + Vl + ¥ -2} +a 2 x(A; 2 )} + 



0(a 3 ) , (33) 



where x(^ 2 ) is a monotonically increasing positive function. Then, we obtain the exact GSE, 
EM and ANVP of the one- dimensional polaron up to the a 3 -order. The results read 

pn ( 3 A 2 /, 63y^-19y^ 29v/6\ 3 n , 4^ 
£ 3 (a) = -«-^=-lJ a -|5 —j a -0(a), 

m* 3 (a) = 1 + l a + (—= -\\ a 2 + 0.0691096281 a 3 + 0(a 4 ) , (34) 



2 \8V2 4 

V2 ~ 2 ) " ~ 2 ^ 16 48"^ 

A refined weak-coupling method based on the three-phonon correction has resulted in 



Ar , , 1 /3 \ 2 7 / 63^-19^ 29 v / 6\ , 4N 
jV 3 (a) = - a 4 (_-2|a 2 --(5 — ^— ) a 3 + 0(a 4 ^ 



S cw (a) = -a - 0.06066a 2 - 0.00844a 3 , 

m cww( a ) = 1 + 0.5 a + 0.19194 a 2 - 0.06912 a 3 , (35) 
Ncww(a) = 0.5 a + 0.12132 a 2 + 0.02954 a 3 . 

Our weak-coupling results fl3"4] ) improve the known data R |7| on the one-dimensional 
polaron. Note, a recent attempt by considering an additional confining dot potential has 
resulted in the binding energy and effective mass diverging in one dimension [|T^]. 

Strong coupling 

In the strong-coupling limit the solution of (PI|) looks like 

D(k) = (k 2 + /j 2 )- 1 , /i = 4a 2 /7r. (36) 

The Gaussian leading-order and the corrected estimates for the GSE, EM and ANVP become 

E (a) = -(l/vr)a 2 , m*(a) = (16/vr 2 ) a 4 , iV (a) = (2/tt) a 2 , 

E 2 (a) = -0.327014a 2 , m* 2 (a) = 1.858065 a 4 , JV 2 (a) = 0.654028 a 2 , (37) 

E 3 (a) = -0.330205 a 2 , m*(a) = 1.966430 a 4 , iV 3 (a) = 0.660410 a 2 . 

Rigorous studies of the strong-coupling regime of the polaron in one dimension have been 
performed, particularly, in [f|, |5]. Within the Pekar adiabatic approximation |7j efficient for 
a — > 00, one arrives at the following variational task with a normalized wave function \l/(x): 

E A {a) = -a 2 Jdx{^(x)A^(x) + 2\^(x)\ 4 } , m* A (a) = a 4 -Jdk k 2 | (^\e ikx \^) \ 2 . (38) 

A modified polynomial placed instead of ^f(x) in @ has resulted in E PS = —0.3330877a 2 , 
m* PS = 2.1254 a 4 and N PS = 0.6661754 a 2 . Fortunately, problem (|38| ) admits explicit ana- 
lytic solutions (obtained first in [|J with a re-scaled coupling) as follows: 

^(x) = V2/(e x + e~ x ), E A (a) = -a 2 /3 , m A (a) = 32 a 4 /15 , N A (a) = 2a 2 /3. (39) 
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We see that our Gaussian (leading-order) approximations E (a), m^{a) and N (a) differ 
from the adiabatic results (^). However, these gaps have been systematically narrowed by 
involving higher-order non-Gaussian contributions in (|37]). 

Intermediate coupling 

For finite a we have solved numerically the constraint equations in (f21~1). The obtained 
numerical results for the Gaussian leading-order contributions E Q (a) and m1j(a) as well as 
for the next-to-leading corrected data are plotted in Fig. 1 in comparison with the Feynman 
estimates ||. In Fig. 2 we have represented N (a) calculated by formula (|^). According 
to our results, the number of virtual phonons increases fast as a grows and hence, the 
three-phonon correlation method |§ should be used carefully for a > 2.5. 



0.01 r 1 0.15 




-0.04 r -0.15 

0346802466 

a a 

Fig. 1 

The ground-state energy and effective mass of the one-dimensional polaron 
normalized to the corresponding Feynman estimates (dotted horizontal lines) 
versus the coupling constant a. Dashed lines represent the leading-order Gaus- 
sian approximations and solid curves correspond to the corrected results taking 
into account the second-order non-Gaussian contributions. 

Note, that the Feynman-type estimation |J represents^ partial case of the Gaussian 
leading-order approximation, when a trial convex function Dp(k) = w/k 2 + (l — w)/(k 2 + {i 2 ) 
is used instead of the exact solution D{k). Therefore, the Feynman estimate is less exact 
than the general Gaussian result, in particular, Eqf(oi) > E (a). 

In conclusion, the ground-state characteristics of the large optical polaron confined in one 
dimension have been calculated by using the generalized Gaussian approximation method. 
Exact analytic solutions for the GSE, EM and ANVP have been obtained up to the order a 3 
in the weak-coupling limit. In the strong coupling limit, our Gaussian leading-order results 
deviate (e.g., by 0.9 percent for the energy) from the adiabatic data obtained within the 
Pekar Ansatz, but the gaps have been systematically narrowed by involving the higher-order 
corrections. For finite coupling, the leading-order Gaussian GSE improves the Feynman vari- 
ational estimate and represents the lowest upper bound available for the one-dimensional 
polaron. Our result for the number of excited phonons limits the validity region of the few- 
phonon approximation methods. The next-to-leading non-Gaussian contributions have been 
calculated to correct the Gaussian approximations for finite a. 



12 3 

a 

Fig. 2 

The Gaussian approximation to the average number of virtual phonons iVo(a) 
versus the coupling constant a. 

Applied to the polaron confined in one dimension, the GGA method provides reliable 
results rather quickly in the entire range of coupling. Therefore, it will be reasonable to 
apply the GGA method to improve the Lieb-Yamazaki type lower-bound estimate (see, e.g. 
||) of the polaron self-energy. 
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